Anisotropy in the Incommensurate Spin Fluctuations of Sr 2 Ru04 
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It has been proposed that S^RuCU exhibits spin triplet superconductivity mediated by ferromag- 
netic fluctuations. So far neutron scattering experiments have failed to detect any clear evidence 
of ferromagnetic spin fluctuations but, instead, this type of experiments has been successful in con- 
firming the existence of incommensurate spin fluctuations near q = (| | 0). For this reason there 
have been many efforts to associate the contributions of such incommensurate fluctuations to the 
mechanism of its superconductivity. Our unpolarized inelastic neutron scattering measurements 
revealed that these incommensurate spin fluctuations possess c-axis anisotropy with an anisotropic 
factor Xc/x'a b °f ~ 2.8. This result is consistent with some theoretical ideas that the incommensu- 
rate spin fluctuations with a c-axis anisotropy can be a origin of p-wave superconductivity of this 
material. 



I. INTRODUCTION 

Sr2Ru04 is the first 2D perovskite oxide material known to exhibit a superconducting transition without containing 
copper 1 . While Sr2Ru04 is isostructural with the high-T c material I^-^Sr^CuO^ its normal state shows Fermi liquid 
behavior and it's superconducting state is not a spin singlet (S — 0) as observed in the conventional s-wave (I = 0) 
superconductors or the d-wave (I = 2)high-T c materials. Its superconducting state is instead a spin triplet (S = 1) 
with (most-likely) p-wave symmetry (1 = 1) (see Ref^ for a recent review). 

/iSR and NMR (Knight shift) measurements have provided experimental evidence of the spin triplet pairing in 
Si'2Ru04. /iSR measurements succeeded in confirming the existence of the spontaneous magnetic field below the su- 
perconducting transition temperature T c ~ 1.5 K, indicating the time- reversal symmetry-breaking in superconducting 
state*. Knight shift measurements for the oxygen site in the Ru02 planes revealed that the spin susceptibility remains 
temperature independent even below T^. 

Since T c of Sr2Ru04 ( ~ 1.5 K) is strongly suppressed by nonmagnetic impurities^,a mechanism other than electron- 
phonon interaction was proposed as the origin of the pairing mechanism of the superconductivity observed in this 
system. From the analogy with the super-fluid state of 3 He and from the fact that the neighbor material SrRuC>3 
is ferromagnetic, it was speculated that S^RuC^ exhibits spin triplet superconductivity mediated by ferromagnetic 
fluctuations. Up to this date, however, there is no clear experimental evidence of ferromagnetic fluctuations in this 
materialSiS. 

The electronic structure of Sr2Ru04 (cJ-electron system) is much simpler than those of other spin triplet 
superconductor a 9 i 10 . This fact has stimulated theoretical efforts on the topics of spin triplet superconductivity and 
the symmetry of superconducting order parameters. 

Mazin and Singh have calculated the electronic band structure of Sr2Ru04 based on the t2 g orbital of the Ru 4+ 
(4d 4 ) and showed that the Fermi surfaces consist of quasi one-dimensional a, (3 planes defined up by the d yz , d zx 
orbital, and of two-dimensional 7 planes defined the d xy orbitalii. These predictions are consistent with the results 
of dHvA experiments^. Furthermore, the theory predicted that sizable nesting effects in the quasi one-dimensional 
bands (a, /3 planes) may cause the enhancement of the spin susceptibility near the incommensurate propagating vector 



2 



q = (i i 0)i£. Such an enhancement was indeed confirmed in dynamical spin susceptibility x"(q ,u;) at q = (0.3 
0.3 0) by inelastic neutron scattering (INS) experimental. These results stimulated discussions about the possibility 
of p-wave superconductivity mediated by such incommensurate spin fluctuations. 

Some theoretical works reported that, if such an incommensurate spin fluctuations possess c-axis anisotropy, the spin 
triplet superconductivity could be stabilized by such fluctuationsi^*iSii&. It is therefore of great importance to stablish 
if there is any observable anisotropy in x"(q ) that can be related to the origin of the spin triplet superconductivity 
in Sr2Ru04. 

Ishida et al. have reported the observation of the anisotropic behavior of the spin susceptibility measured by the 
NMR techniqueii. In NMR measurements, one can observe the q-integrated spin susceptibility, ^„ x ^ |a>->0- 
Judging from the similarities with the INS data reported by Sidis et al. 6 , Ishida et al attributed the temperature 
dependent part of the q-integrated spin susceptibilities (observed by the NMR) to the spin susceptibility at q , and 
reported that x" (q , uj) has a c-axis anisotoropy with an anisotropic factor x" c /x"a,h, ( x"/c,out(qo7 LJ )/x"/c,m(qo) w ) 
in their notation) of ~ 3. 

In order to ascertain the anisotropic nature of the incommensurate spin fluctuations in Sr2Ru04, however, it is 
necessary to measure its q dependent spin susceptibility x"(q, u>) using INS. We have performed such measurement and 
found that the dynamical spin susceptibility of this system at q indeed exhibits c-axis anisotropy with an anisotropic 
factor of ~ 2.8. This value is consistent with the anisotropic factor estimated from the NMR measurements^ . Our 
conclusion is different from those of the recent reports by Servant et al. and Braden et alm&, this will be discussed at 
the end of this paper. 

II. EXPERIMENT 
A. Sample preparation and Experimental setup 

For this neutron scattering study, we grew large single crystals of Sr2Ru04 by the floating zone method. The 
crystals were cut into smaller cylindrical pieces (4 mm</> in diameter and 30 mm in length). We performed resistivity 
measurements on these crystals using a Quantum Design PPMS instrument equipped with a 3 He option. These 
measurements revealed that T c (onset) of all the samples lies between 1.4 ~ 1.6 K. 

The unpolarized INS experiments were performed using the triple axis spectrometer GPTAS installed at the JRR- 
3M reactor at the Japan Atomic Energy Research Institute (JAERI) in Tokai, Japan. Neutrons with a fixed final 
momentum of kf = 3.83 A -1 and a combination of horizontal collimations of 40'-80'-40'-80' (FWHM from the 
monochromator to the detector) were utilized. A pyrolytic graphite (PG) filter was placed after the sample position 
to eliminate higher order wavelength contaminations. Three sets of crystals were prepared in order to probe three 
different scattering planes, (h k 0), (h h I) and (0.7ft- 0.3h I). The total volume of each of these sets was ~ 3 cm 3 —. 
The crystals were sealed in aluminum cans (filled with He exchange gas to ensure a uniform temperature) that were 
attached to the cold head of a closed-cycle He gas refrigerator. Throughout this paper the scattering vector Q = (Qh 
Qk Qi) is indexed in reduced lattice units with tetragonal reciprocal lattice of a* = b* = 1.63 A -1 and c* = 0.49 A" 1 , 
respectively. 

B. Magnetic neutron scattering 

In this section we describe the method that we used to measure the anisotropic factor of the spin susceptibility 
x"(q, uj). In a magnetic neutron scattering experiments, the scattering intensity I is given by 

I oc f Q x G, (1) 

where /q is the magnetic form factor, which is the Q component of the fourier transform of the distribution of 
unpaired electrons that contribute to the magnetism in the system. If the electronic distribution is isotropic, /q 
shows a monotonic decrease with the absolute value of the scattering vector Q, Q as demonstrated in Fig. 1(a). The 
quantity G in Eq.(l) is an orientation factor related to the fact that neutrons are scattered only by the magnetic 
components perpendicular to the scattering vector Q. In the present study, we assumed that the spin susceptibility 
within the Ru02 planes in tetragonal Sr2Ru04 is isotropic (x" = X& = x" b)- The orientation factor G is then given 
by 



G(9) = (1 + sin 2 6) xl', fc (q, w) + cos 2 6 x£(q, to), 



(2) 
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where 6 is the angle between the scattering vector Q and the ab plane, which changes through the Qi component of 
the scattering vector. 

In Fig. 1(b) we show a calculation of the Q dependence of G(8) for Q=(0.3 0.3 Q{) for different anisotropic factors. 
G(0) is constant when the susceptibility is isotropic (x" b = x") but shows different Q dependence with anisotropy 
(Xab Xc)- Namely, G(9) increases (decreases) with Q when Xc/Xab < 1(> ■*-)■ 

The Q dependence of the total intensity I given by Eq.(l) is shown for (0.3 0.3 Qi), (0.7 0.3 Qi) and (0.7 0.7 Qi) in 
Fig. 1(c). If x"(q, u) is isotropic, the intensity / is scaled only by fg but it decreases slower (faster) than fg in the 
presence of anisotropy Xc/x'a b < 1(> !)■ This illustrates the fact that the anisotropic nature of the spin fluctuations 
can be directly determined by the comparison of the Q dependence of the intensity / and fg. We would like to stress 
that knowing the magnetic form factor accurately is the key to the accurate determination of the spin susceptibility 
anisotropy factor. Unfortunately the only magnetic form factor that can be found in the literature for the Ruthenium 
is that for Ru + (/q(Ru + ))SS. The Ru ions in Sr 2 RuC>4 are not Ru + but Ru 4+ (nominally) . Furthermore, Sr 2 RuC>4 
is not an insulator but an itinerant electron system, and the use of /q(Ru + ) to characterize the magnetic response 
of Sr2RuC>4 is clearly inadequate. For this reason we decided to determine the magnetic form factor for Sr2Ru04 
(/q(Sr2Ru04)) experimentally. 

C. Determination of the magnetic form factor of Sr2RuC>4 

To determine the magnetic form factor of S^RuCU, /g(Sr2Ru04), we measured the Q dependence of spin suscep- 
tibility at several Q positions with q =(0.3 0.3 0) in the (h k 0) plane (8 = 0). The Q dependence of the observed 
intensities is shown in Fig. 2, the filled and open symbols indicate our data and those reported data by Sidis et al&, 
respectively. Note that, throughout the present paper all the quoted intensities have been corrected for resolution- 
volume effects, and that all the quoted experimental errors correspond to 2 x a in order to reflect the ambiguities of 
the scattering technique. 

Note that (0.7 0.3 0) and (1.3 0.3 0) are not q positions from the reciprocal zone center, T point, but those from 
the Z point (ex. (1 0)). These data can be treated equally with other data, because the spin susceptibility at q 
shows a strong two dimensionality and a rod type scattering along the c*-axis so that one can observe the signal even 
on the (h k 0) zone. 

Our first observation of Fig. 2, is that it is clear that the data do not scale with /q(Ru + ), and decrease faster 
than it. This behavior is consistent with the fact that S^RuC^ is an itinerant electron system where mobile electrons 
distribute wider in real space than localized electron system and strongly indicates that one can not use /q(Ru + ) to 
evaluate anisotropic factor of spin susceptibility of Sr2Ru04. 

The full line in this figure is /g(Sr 2 Ru04) determined in the present study, for this purpose we fitted the observed 
intensities to the expression 

/ Q (Sr 2 Ru0 4 ) = A exp [B (Q/4vr) 2 ] + C. (3) 

Here, we assumed that the /q in the (h k 0) plane is isotropic, so that the Q dependence of the /q is described as a 
single Q functional. 

This form factor was used to evaluate the anisotropic factor of the incommensurate spin fluctuations. Note that 
conductivity and coherence length of Sr 2 Ru04 show anisotropic behavior (a a ,b > &c and £_ a ,b > £c) 2 - Such results 
indicate that a distribution of unpaired electrons along the c-axis is confined and then the decrease of /q with Qi 
must be slower than that for Qh or Qk- It should be stressed here that we assumed an isotropic form factor /q in all 
directions in the present study which causes an underestimation of the c-axis anisotropy. 

III. EXPERIMENTAL RESULTS 

A. Qi dependence of intensity 

In order to study the Qi (8) dependence of the intensity, we performed several series of constant-i? scans along 
(0.3 0.3 Qi), (0.7 0.7 Qi), (0.7 0.3 Qi) and (1.3 0.7 Qi), and found that, because of low intensity, it is difficult to 
get accurate Qi dependence at Q positions farther than (0.7 0.7 Qi). For this reason we report only the results at 
the (0.3 0.3 Qi) and (0.7 0.3 Qi) scans. Furthermore, to collect reliable data, one needs to select a clear window of 
energy where any spurious peaks including phonon scattering do not appear. For the constant-E' scans at (0.3 0.3 
Qi) and (0.7 0.3 Qi), neutron transfer energies were selected to be 4 meV and 8 meV, respectively, by measurements 
of energy dependence of intensity at (0.3 0.3 0) and (0.7 0.3 0) with energy transfer between ~ —2 < E <~ 20 meV. 
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The energy dependence of intensity at (0.3 0.3 0) is shown in an inset of Fig. 3(a). The result clearly shows that the 
intensity at E = 4 mcV is affected by neither incoherent nor forward scattering. 

Qi dependence of integrated intensity at (0.3 0.3 Qi) and (0.7 0.3 Qi) are depicted in Fig. 3(a), the integrated 
intensities were calculated as the product of intensities at (0.3 0.3 Qi), (0.7 0.3 Qi) and the width determined by 
constant- E scans along the (h k 0) direction. The obtained widths at (h h Qi) and (0.7h 0.3h Qi) were almost constant 
with Qi, and we used their averaged values, cf. 0.188 and 0.184 A -1 (in FWHM), respectively. In addition, intensities 
at (0.5 0.5 Qi) and (0.7±0.1 0.3±0.04 Qi), which are almost constant with Qi, were subtracted as background for 
calculations of peak intensities at (0.3 0.3 Qi) and (0.7 0.3 Qi), respectively. Finally the data at (0.7 0.3 Qi) with 8 
meV were scaled with the data at (0.3 0.3 Qi) with 4 meV by detailed measurements of energy dependence of signals. 

The Qi dependence of the integrated intensity (0.3 0.3 Qi) and (0.7 0.3 Qi) in the Fig. 3(a) shows a very broad 
peak centered at Qi = , indicating the strong two dimensionality of the spin fluctuations. This result is consistent 
with the one reported by Servant et ali, and allows us to neglect the magnetic correlations along the c-axis. Thus 
we treat data sets with different Qi independently. 

B. Determination of the anisotropic factor 

Fig. 3(b) shows the Q dependence of the intensities for (0.3 0.3 Qi) and (0.7 0.3 Qi), the full line is the magnetic 
form factor /^(S^RuC^) that we measured as indicated above. This figure clearly shows that the intensities for both 
(0.3 0.3 Qi) and (0.7 0.3 Qi) decrease faster than /|(Sr 2 Ru0 4 ) with increasing Q. Such Q dependence correspond to 
the case with x" b < Xc as demonstrated in Fig. 1(c). 

To evaluate the anisotropic factor, Q dependence data for (0.3 0.3 Qi) was fitted to Eq. (1)(2), the data for (0.7 
0.3 Qi) was not used in this fit because of the poor statistics. From the fitting, we calculated the anisotropic factor 
of the spin susceptibility, Xc/x" b ~ 2. 8 ± 0.7. Note that to evaluate an error in the determination of the anisotropic 
factor we took into account the error of the magnetic form factor. Furthermore, as explained in the previous section, 
we assumed an isotropic form factor and such assumption may cause the underestimation of the c-axis anisotropy. 
These results let us conclude that the incommensurate antiferromagnetic fluctuations observed at q = (0.3 0.3 0) 
exhibit c-axis anisotropy. 

IV. DISCUSSION 
A. Static and dynamical spin susceptibility observed in Sr2Ru04 

The magnetic properties in the normal-state of Sr2Ru04 reported so far are (A) slightly anisotropic uniform 
susceptibility at q = 2 -, (B) anisotropic spin fluctuations with anisotropic factor of Xclx'a & ~ 3 at somewhere in q 
reported by NMRii, (C) spin fluctuations observed at incommensurate q of (0.3 0.3 0) observed by INS 6 . 

The uniform susceptibility of (A) is explained by the Pauli paramagnetism of the conduction electrons in a two- 
dimensional 7 band, and the origin of a slight anisotropy of it (xc/Xa.b ~ 1-1) is attributed to the orbital Van Vleck 
contribution, which is affected by fields parallel to the c-axis due to the 1-dimensional a- and /3-bands 23 . 

On the other hand, the anisotropy in (B) can not be associated with that in (A), because the anisotropic factor and 
energy scale of each spin susceptibilities are quite different. Judging from the similarity in temperature dependence 
of spin fluctuations in (B) and (C), Ishida et al. speculated that the anisotropic behavior observed in the NMR 
measurements has a close relation with spin fluctuations observed at (0.3 0.3 0)^£. Supporting this, our result 
clearly revealed that the incommensurate spin fluctuation has anisotropy with an anisotropic factor of Xc/x'a b ~ ^.8. 
The anisotropic factor reported by NMR measurement is ~ 3, which is in good agreement with the present result. 
These results let us conclude that anisotropic behavior observed by the NMR measurements is associated with spin 
fluctuations at incommensurate qo vector of (0.3 0.3 0). 

B. The origin of the anisotropic behavior 

A short comment about the origin of the anisotropy of the incommensurate spin fluctuations observed at q = (0.3 
0.3 0). Theoretical calculations within the random-phase approximation using a three-band Hubbard Hamiltonian 
predict that spin-orbit coupling plays an important role and that, due to strong coupling, the out-of-plane component 
of the spin susceptibility is about two times larger than the in-plane one at low temperature 2 ^. The calculated 
anisotropy and our result are quantitatively consistent. Magnetic property of (A), (B), and (C) deeply connect with 
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the orbital of d-electrons in R11O2 planes. These facts strongly indicate that the spin-orbit interaction is important 
to discuss the magnetic properties of this system. 

C. Relation between incommensurate spin fluctuations and the superconducting mechanism 

As described in the introduction, some theoretical groups reported that incommensurate spin fluctuations with a 
c-axis anisotropy, x'c > Xabi may introduce a spin triplet superconducting state and that the d vector turns to a 
direction of larger antiferromagnetic fluctuationsi^S^. Our results show that the incommensurate spin fluctuations 
observed in Sr2Ru04 satisfy this requirement, namely x'c > Xo6> which makes a direction of d vector to be parallel to 
the c-axis consistent with the experimental observations^^. 

Then the question here is whether these spin fluctuations are really driving forces of the superconductivity of this 
material or not. Basically, the superconductivity of Sr2Ru04 is believed to originate in the quasi two-dimensional 
7 main band. On the other hand, the incommensurate antiferromagnetic fluctuations of Sr 2 Ru04 is caused due 
to the nesting property of the one-dimensional a and (3 bands. Furthermore, in the Sr2Rui_ x Ti x 04 (in which 
superconductivity is quickly suppressed and the antiferromagnetic fluctuations observed at x = develop into a static 
order with increasing x) , the x dependence of T c seems to be explained only by a doping effect and no enhancement of T c 
by the spin fluctuations was observed2£i2L2£. These results imply that the incommensurate spin fluctuations may not 
contribute to its superconducting mechanism^. In order to further clarify the mechanism of the superconductivity in 
Sr 2 Ru04, especially of relations between spin triplet superconductivity and antiferromagnetic fluctuations, information 
of x"(q 0) cj) behavior below T c would be of great help. 

D. Discrepancies with unpolarized INS results by other groups 

In the present study, we measured data at (0.3 0.3 Qi) and (0.7 0.3 Qi) including (0.3 0.3 0) (Q ~ 0.70 A -1 ) 
and estimated the anisotropic factor x'c/x" b t° be ~ 2.8 by evaluating a difference between the / and /g(Sr2Ru04). 
These results, however, are at odds with other unpolarized neutron scattering data reported by Servant et ali and 
by Braden et alA. We attribute these discrepancies to (a) the narrower Q range in these groups' measurements, and 
(b) a lack of a determination of /Q(Sr2Ru04) by the other groups. For example, Servant et al. measured data at 

(0.3 0.3 Qi) with only larger Q part (Q > 0.80 A -1 ) and (0.7 0.7 Qi), and concluded an isotropic behavior of spin 
fluctuations based on the fact that a small number of data points observed along (0.3 0.3 Qi) and (0.7 0.7 Qi) scaled 
at a very narrow Q range near Q of ~ 1.6 A -1 . One can see in Fig. 3(b) of our paper, that the accuracy of the data 
in that Q range (near ~ 1.6 A" 1 ) is not very good. We also observed data at (0.7 0.7 Qi) and found that the data 
scaled with those along (0.3 0.3 Qi) in this Q range within the huge error bars. Furthermore, they did not get the 
proper magnetic form factor for Sr2Ru04, this prevented them from making a reliable comparison of their data with 
the magnetic form factor in the small Q region. On the other hand, Braden et al. observed Q-dependence of I at 
(0.3 0.3 Qi) within a very narrow Q range of 1.2 ~ 2.5 A -1 and showed that the data decreases slower than that of 
/q(Ru + ) 8 . This behavior is clearly opposite to the our data shown in this paper and to the data of Servant et al, 7 . 
At this time we do not understand the source of this discrepancy. 

Recently, neutron polarization analysis experiments have been performed on Sr2Ru04 by two independent groups. 
These groups succeeded in confirming the c-axis anisotropy with anisotropic factor of 2^2.52^ and 2.0 ± 0.42S), 
respectively, being consistent with our unpolarized neutron results presented here. It can be argued that the best 
way to perform this type of measurements is using the polarized neutron scattering technique because this technique 
allows the separation of magnetic components to the scattering from any other non-magnetic components including 
phonon and spurious peaks. We would like to stress here, however, that this is not the only reliable way to measure 
magnetic components. It is true that the unpolarized neutron technique intrinsically more ambiguous when it comes 
to measure magnetic components. But being conscious of this fact, we paid the greatest care to reduce such errors 
and we made many consistency checks with different scattering zones and even checked background from the cryostat 
and judiciously chose the best conditions for the experiment. Every experimental result presented in our paper has 
been examined with great caution and our results are reliable. 

V. SUMMARY 

We have performed unpolarized inelastic neutron scattering measurements on Sr2Ru04 to probe the anisotropic 
behavior of the spin susceptibility observed at the incommensurate wave vector of q = (0.3 0.3 0). Our measurements 
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indeed support that the susceptibility exhibits a c-axis anisotropy i.e. x'c/x" b ~ 2.8 ± 0.7. This anisotropy ratio is in 
good agreement with the result obtained by the NMR measurements (~ 3)A£. 
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FIG. 1: (a) A schematical Q dependence of the square of magnetic form factor, fg. (b) Q dependence of calculated orientation 
factors G(6) at Q = (0.3 0.3 Qi) for different anisotropic factors. Depending on a ratio, x" b '■ Xc, G(8) shows different Q 
dependence, (c) Q dependence of calculated intensities I (= fg x G(8)) at (0.3 0.3 Qi), (0.7 0.3 Qi) and (0.7 0.7 Qi) with 
different anisotropic factors. 




FIG. 2: Q dependence of intensities observed at Q with q =(0.3 0.3 0) in the (h k 0) plane (9 = 0). Filled symbols are the 
present results, in which circle and diamond symbols correspond to the data taken with different sample set with (h k 0) and (h 
h I) with I — 0, respectively. Open symbols are taken from Refp2> reported by Sidis et ai. Taking into accounts ambiguities of 
scattering experiments, we conservatively adopt 2 x a error bars. A broken and a full lines correspond to the square of magnetic 
form factor of Ru + , /q(Ru + ), and that of S^RuCU, /^(S^RuCU), respectively. The latter was determined and parameterized 
with Eq. (3) in the present study (See the main text). 
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FIG. 3: (a) Qi dependence of integrated intensity / at (0.3 0.3 Qi) and (0.7 0.3 Qi) and at 4 K. Some data points are missing 
because of uncertainty due to existence of spurious peaks, (inset) : Energy dependence of scattering intensity at (0.3 0.3 0). 
(b) Q dependence of resolution corrected intensity at (0.3 0.3 Qi) and (0.7 0.3 Qi). Averaged values between (Qh Qk \Qi\) and 
(Qh Qk — \Qi\) are plotted. The full line is the square of magnetic form factor of S^RuCU, /^(S^RuOzi), determined in the 
present study (See Fig. 2). Dashed line on (0.3 0.3 Qi) is a fitting curve to Eq. (1)(2) with x'c/x'a,b — 2.8 and a broken line 
on (0.7 0.3 Qi) is a simulation line calculated with the parameters evaluated by data at (0.3 0.3 Qi). In both plots, 2 x a was 
conservatively adopted as error bars. 



